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ABSTRACT 


This  report  investigates  the  use  of  energy  bounds  to  estimate:  the  amount  of  energy  that 
can  be  transmitted  through  an  aperture  when  exposed  to  an  EM  pulse;  the  fields  inside  a  cavity 
with  one  or  more  apertures;  and  the  coupling  of  energy  to  a  wire  inside  a  cavity  with  apertures. 
New  bounds  have  been  determined  for  the  energy  that  can  be  coupled  from  a  damped-sine  wave, 
a  rectangular  RF  pulse  and  an  inverse-exponential  pulse. 

It  is  shown  that  the  use  of  energy  bounds  is  a  useful  approach  to  the  estimation  of  interior 
fields  and  coupling  to  systems  inside  a  structure  with  apertures.  This  approach  has  potential  for 
simplifying  the  hardening  of  electronic  systems. 


RESUME 

Ce  rapport  examine  I’utilisation  des  limites  superieures  d’energie  pour  determiner:  la 
quantite  d’energie  tramsmise  par  des  ouvertures  lorsqu’elles  sont  exposes  a  des  impulsions 
electromagnetiques;  les  champs  electromagnetiques  a  I’interieur  d’une  cavite  avec  des 
ouvertures;  et  la  couplage  a  un  cable  a  I’interieur  d’une  cavite  avec  des  ouvertures.  Les  limites 
superieures  d’energie  ont  ete  developees  pour  des  impulsions  sinusoidale  amortie,  des  signaux 
RF  rectangulaire,  et  des  signaux  exponentiel. 

On  a  demontre  que  I’utilisation  des  limites  superieure  est  une  approche  efficace  pour 
determiner  les  champs  interieur  et  la  couplage  a  des  systemes  a  I’interieur  d’une  structure.  Cet 
approche  peut  etre  efficace  pour  le  durissement  electromagnetique  des  systemes  electronique. 
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EXECUTIVE  SUMMARY 


The  electromagnetic  (EM)  vulnerability  of  CF  vehicles,  weapons  systems, 
and  communications  systems  is  of  concern.  Potential  threats  include  natural  sources,  such  as 
lightning  and  electrostatic  discharge,  as  well  as  a  variety  of  man-made  sources  that  includes  the 
"normal"  operational  EM  environment  produced  by  high  power  radars,  EW  transmitters,  and 
communications  systems.  At  the  same  time,  technological  advances  such  as  the  use  of  a 
composite  materials  for  construction;  and  increased  miniaturization  of  electronics  are 
contributing  to  increased  EM  susceptibility  of  many  systems. 

Analysis  of  the  effects  of  external  fields  on  military  electronic  systems  is  an  important 
part  of  the  hardening  process.  For  large  systems,  such  as  a  ship  or  aircraft,  extensive  computing 
facilities  are  needed  to  calculate  the  interior  fields  inside  even  a  simplified  model.  In  this  report, 
an  alterative  approach  based  on  the  use  of  energy  bounds  is  examined.  It  is  shown  that  this 
method  has  potential  for  simplifying  the  EM  hardening  of  military  systems. 
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1.  BACKGROUND 


A  knowledge  of  the  interaction  and  coupling  of  electromagnetic  fields  with  large,  closed 
structures  such  as  aircraft,  ships,  missiles  and  buildings  is  critical  for  a  designer  to  ensure  that 
installed  electronic  equipment  is  compatible  with  the  external  electromagnetic  environment.  The 
external  interaction  problem  usually  involves  the  calculation  of  the  surface  currents  induced  on 
a  simplified  wire  or  patch  model  of  the  structure  and  subsequent  calculation  of  the  interior  fields 
[1] .  When  the  structure  become  electrically  large,  the  calculation  of  the  induced  surface  currents 
and  fields  requires  extensive  computer  resources  even  when  the  model  is  simplified  to  the 
maximum  extent  and  maximum  use  of  symmetry  is  used.  For  example,  calculations  of  the 
penetration  of  fields  into  the  bridge  of  the  CPF  [1]  took  about  18  hours  on  a  VAX  6000  at  each 
frequency  and  was  limited  to  a  maximum  frequency  of  about  200  MHz. 

Because  of  the  requirement  for  extensive  computing  facilities  to  calculate  the  fields  inside 
even  a  simplified  model  of  a  large  structure,  it  is  usually  not  possible  to  include  interior 
equipment  in  the  calculation.  When  this  is  not  possible,  the  coupling  of  fields  to  electronic 
systems  is  normally  carried  out  in  two  stages.  The. first  step  is  to  calculate  the  interior  fields 
using  a  simplified  model  that  neglects  the  interior  components.  The  second  step  is  then  to 
calculate  the  coupling  of  these  fields  to  the  electronic  systems  of  interest  and  to  estimate  system 
susceptibilities. 

The  above  approach  has  the  following  difficulties: 

-  if  there  is  strong  coupling  to  the  interior  fields,  the  Q  of  the  structure  will  be 
lowered  and  the  fields  calculated  will  be  overestimated; 

-  the  coupling  calculated  using  this  two-stage  approach  can  violate  energy 
conservation  (ie.  the  calculated  coupled  energy  is  more  than  that  which  can  be 
transmitted  through  the  apertures  in  the  structure). 

-  neglect  of  minor  features  often  gives  models  a  higher  symmetry  than  what  exists 
in  reality  in  order  to  reduce  computation  time.  Because  these  models  are 
expected  to  have  higher  Q -factors,  the  interior  fields  may  be  overestimated. 
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Because  of  the  difficulties  outlined  in  the  computational  approach  described  above,  an 
alternative  approach  based  on  the  use  of  energy  bounds  has  been  proposed  by  Lee  [2],  Wame 
and  Chen  [3,4]  and  others  [5].  In  this  approach,  one  attempts  to  determine  the  maximum 
energy/power  that  can  couple  to  an  electronic  system  inside  a  structure  with  apertures.  This 
approach  recognises  that  it  is  virtually  impossible  to  model  all  of  the  detail  of  any  real  world 
system  and  attempts  to  find  a  solution  (in  terms  of  an  upper  bound)  that  is  independent  of  these 
system  details. 

In  the  remainder  of  this  report,  the  use  of  bounds  is  explored  as  a  method  of  estimating 
interior  fields  and  coupling  to  systems  inside  of  a  structure  with  apertures. 
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2.  BOUNDS  ON  ENERGY  TRANSMITTED  THROUGH  AN  APERTURE 


2.1  APERTURE  THEORY 


Consider  the  penetration  of  the  electromagnetic  field  through  an  aperture  in  an  infinite 
planar  screen,  assumed  to  be  perfectly  conducting  and  vanishingly  thin  ,  as  shown  in  Figure  1 . 
Previous  work  has  shown  that  this  problem  is  most  conveniently  treated  using  the  equivalence 
principle  [6].  Using  this  principle,  the  penetration  of  the  fields  through  the  aperture  is 
determined  by  considering  the  screen  to  be  continuous  and  by  placing  appropriate  magnetic 
surface  currents,  M^,  over  the  aperture  (Figure  2).  The  magnetic  field  {H~)  on  the  left 
(illuminated)  side  can  be  written  [7] 

H~ir)  =H‘^-—[k'^Fir)  +V(V*F(f ))].  (1) 

In  this  expression,  F(f )  is  the  electric  vector  potential  given  by 

F(f )  =^// 2Miif)G(rf')ds' 


where 


(3) 


and 


G(r/')  = 


g  -ik\r-r'\ 


(4) 


In  these  equations,  is  the  short  circuit  magnetic  field  [7],  k  is  the  propagation 
constant  and  o  is  the  angular  frequency. 
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(5) 


On  the  right  hand  side  of  the  screen  ,  the  magnetic  field,  ,  is  given  by, 

H\n  =  •^[^W)+V(V-F(f))]. 

The  electric  fields  in  the  two  half  spaces  are  given  by  the  relationships, 

E~ir)  =  E‘^  -  -VxF(f)  (6) 

€ 

and 

E\r)  =  -VxF(f)  (7) 

€ 


where  F“  is  the  short  circuit  electric  field  [7]. 

When  the  aperture  is  open,  the  transverse  components  of  the  magnetic  fields  are 
continuous  and  the  following  relationship  can  be  written  [7] 

lim(^’(f)x  «,  )  =  lim(^*(f)x  u  )  (8) 


which  from  equations  (1)  and  (5)  leads  to  the  relationship 

i  =  — [V,V*F(f)+;t2F(f)]  X  a  .  (9) 

it*  '  '  ^ 


These  equations  show  that,  by  using  the  incident  magnetic  field,  H\  the 
equivalent  magnetic  currents,  ,  can  be  found  by  solving  (9)  and  the  EM  fields  in  both  half 
spaces  can  be  determined. 

2.2  SMALL  APERTURE  THEORY 

The  electromagnetic  field  that  penetrates  through  an  electrically  small  aperture  can  be 
approximated  [8]  by  the  radiation  from  a  magnetic  dipole  moment,  ,  and  an  electric  dipole 
moment,  p^,  both  located  at  the  centre  of  the  aperture.  These  dipole  moments  are  related  [8] 
to  the  magnetic  surface  current,  M^,  by  the  expressions 

Pe  =  (10) 

^  A 
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and 


?.  =  — /Me?')*'-  (11) 

0)|iy 

The  above  relationships  show  [9]  that  the  magnetic  and  electric  dipole  moment  are  related 
to  the  irrotational  component  and  the  solenoidal  component  of  respectively. 

The  equivalent  dipole  moments,  and  ,  are  proportional  to  the  incident  field 
strengths  and  and  electric  and  magnetic  polarizabilities,  and  ,  can  be  defined  by 

Pe  =  “z 


and 


Pm  =  • 

^  m  m  t 


(13) 


In  these  expressions,  is  the  normal  component  of  the  incident  electric  field  and^/ 
is  the  tangential  component  of  the  incident  magnetic  field. 


For  a  small,  circular  aperture,  the  electric  and  magnetic  polarizabilities,  and  a^,  are 
scalar  and  are  given  [10]  by 


4  3 

a,  =  — 

*  3 


a. 


8  3 

3 


(14) 


where  a  is  the  radius  of  the  aperture.  Analytical  expressions  have  also  been  given  [11]  for 
small,  elliptical  apertures  and  numerical  results  have  been  given  by  De  Meulenaere  and  van 
Bladel  [12]  for  a  variety  of  aperture  shapes. 
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2.3  APERTURE  CROSS-SECTIONS 


The  total  energy  absorbed  and  scattered  by  an  aperture  is  related  to  the  total  aperture 
cross-section,  o^,  by  the  expression 


W. 


tot 


^  J_  f  |E(m)|" 

L  >io 


Oj(Ci))  d(i> 


(15) 


where  tIq  is  the  impedance  of  free  space. 

It  has  been  shown  by  Bom  and  Wolf  [13]  that  the  total  cross-section,  o^((«>),  is  related 
to  the  electric  vector  potential,  F,  by 

a  (o)  =  ^Re(u,’F)  (16) 

*  »  9  * 


where  Re  denotes  the  real  part  of  the  quantity  in  the  brackets  and  is  the  unit  vector  in  the 
direction  of  propagation  of  the  incident  field.  In  the  far  field,  the  scattered  field,  E^,  can  be 
expressed  in  terms  of  F  by  the  relationship 


Akr 


iF- 


kr 


(17) 


If  the  scattering  amplitude,  5(ci>) ,  is  defined  as 


5(0))  = 


(18) 


then 


Oj(o))  =  4nc^Re(S((j3))  . 


(19) 


In  terms  of  the  scattering  amplitude,  the  total  energy,  can  be  written  as 


K,  =  — /|E(M)|*Ee(S(«))<iw  . 
Bo 


(20) 
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Because  the  scattering  process  is  causal,  it  follows  [14]  that  S(o>)  must  be  free  of 
singularities  in  the  Im(o>)>0  half  space.  In  addition,  causality  also  implies  that  the  real  and 
imaginary  parts  form  a  Hilbert  pair  [15].  Use  of  this  property  gives  [15] 


lm(S(cd))  =  - 


2(0 

n 


fRe(S(0) 

J  t)^  - 


dC  =  - 


to  f  o(C) 
In^c^J  C'  -  (0^ 


(21) 


Taking  the  low  frequency  limit,  o>-0,  leads  [15]  to 

OO 

jo^(X)dX  =  .  (22) 

0 


The  low  frequency  limit  of  (22)  can  be  evaluated  [3]  in  terms  of  the  static  electric  and 
magnetic  dipole  moments  of  the  aperture  and  it  becomes 

OO 

fG,(X)dX  =  2Tz^(a^-aJ  .  (23) 

0 


where  and  are  principal  components  of  the  electric  and  magnetic  polarizability  tensors 
of  the  aperture. 


General  analytical  expressions  for  the  cross-section  of  even  the  simplest  aperture  shapes 
are  not  available.  At  low  frequencies,  the  cross-section  is  of  the  general  form  o,(o))  ~  and 
is  also  directly  related  to  the  electric  and  magnetic  polarizabilities  of  the  aperture.  The  following 
analytical  expressions  for  the  cross-section  of  a  small,  circular  aperture  have  been  given  by 
Bethe  [16]  and  Bouwkamp  [17]. 


cos^0‘ 

27tz 


U 


64 

27jt 


+  -sin^6' 

4 


(24) 


where  6'  is  the  angle  of  incidence  relative  to  a  vector  normal  to  the  aperture. 
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At  high  frequencies,  the  aperture  cross-section  approaches  the  geometric  area  of  the 
aperture  and  is  independent  of  frequency  and  polarization,  thus 

o,  “  i4cos0'  .  (25) 


The  shape  of  the  cross-section  when  the  wavelength  of  the  incident  wave  is  comparable 
to  the  size  of  the  aperture  is  not  well  understood.  Crevier  and  Auton  [5]  proposed  the  following 
general  expression  for  the  aperture  cross-section 


a,(o)  = 


i4 


(0>-0)(P^  +  Oo 


(w+Qq)^  +  00 


(26) 


These  authors  used  this  expression  to  derive  a  bound  on  the  energy  transmitted  through 
an  aperture  with  an  incident  field  in  the  form  of  a  "one  minus  an  exponential"  pulse. 


The  following,  approximate,  piece-wise  linear  expression  has  been  given  by  Lee  and 
Yang  [18]  for  the  cross-section  of  a  narrow  slot. 


0,(f)  =  Ifvif,)' 

-  iXiff 


for  fif  = 
for  M,  . 


c 


(27) 


The  following  section  shows  how  other  workers  (especially  Wame  and  Chen)  have  used 
the  above  methodology  to  define  a  bound  on  the  total  amount  of  energy  that  can  be  transmitted 
through  an  aperture. 

2.4  BOUND  ON  ENERGY  TRANSMITTED  THROUGH  AN  APERTURE 

2.4.1  A  Step-Function  Bound  on  EMP  Coupling  131 

The  step  function  can  be  considered  as  the  limiting  case  of  a  double  exponential  pulse 
where  the  rise  time  approaches  zero  (p-<»)  and  the  fall  time  approaches  infinity  (a-*0).  In  this 
case  the  spectrum  of  the  incident  field,  E(<>>) ,  can  be  written 

E(o))  =  — (28) 
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Substitution  of  the  above  expression  for  £(«)  into  equation  (15)  yields  the  following 
expression 

oo 

271  0 

From  (23)  the  following  upper  bound  is  obtained 

2.4.2  A  Bound  on  EMP  Coupling  for  a  Double  Exponential  Pulse  141 

It  is  clear  that  use  of  the  bound  based  on  a  step  function  will  overestimate  the  energy 
transmitted  through  an  aperture  for  a  double  exponential  pulse  with  a  finite  rise  and  fall  time. 
In  a  recent  report,  Crevier  and  Auton  [5]  investigated  the  effect  of  finite  rise  time  on  the  energy 
that  is  radiated  through  an  aperture.  To  do  this  they  assumed  the  model  for  the  aperture 
scattering  amplitude  presented  in  equation  (26),  in  conjunction  with  a  "one-minus-an- 
exponential"  waveform,  to  determine  a  bound  on  the  energy  radiated  through  an  aperture.  The 
form  of  the  scattering  amplitude  chosen  by  Crevier  and  Auton  does  not  satisfy  the  causality 
requirement  that  it  be  free  of  singularities  in  the  Im(co)>0  half  space  however,  in  spite  of  this 
apparent  deficiency,  the  results  obtained  by  Crevier  and  Auton  are  consistent  with  the  results 
obtained  by  Wame  and  Chen  described  below. 

In  a  recent  paper,  Warne  and  Chen  [4]  have  developed  a  bound  on  EMP  coupling  for  a 
double  exponential  pulse  using  a  more  rigorous  treatment  based  on  the  methodology  described 
in  section  2.3.  A  summary  of  this  method  is  given  below. 

The  following  double  exponential  form  (Bell  curve)  is  often  used  as  an  approximation  to 
the  EMP  waveform  because  of  its  simplicity  in  spite  of  the  fact  that  the  first  derivative  is 
discontinuous. 

E(t)  =  Eo(e-“^  -  (31) 
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The  spectrum  of  this  pulse  is  given  by 


(32) 


£(a>)  =  f 

—00 


£o(P  -  «) 


The  total  energy  scattered  by  and  transmitted  through  the  aperture  is  then  given  by 

=  —Eic\^-afRe 
^0 


(33) 


Wame  and  Chen  [4]  have  shown  that,  by  using  contour  integration,  the  above  expression 
reduces  to 


W. 


tot 


2%c^ 


Eq  (p-g) 

ho  (P+«) 


Re 


S(id) 
.  a 


5(iP) 

P  J’ 


(34) 


2.4.3  A  Bound  on  Coupling  for  a  Damped-Sine  Pulse: 

Another  waveform  of  practical  interest,  and  one  that  has  not  been  treated  previously,  is 
the  damped  sine  wave  described  by  the  following  equation 

E(t)  =  Eq  sinco^r  e  u(t)  (33) 

The  spectrum  of  this  pulse  is  given  by 

£(0))  =  fef)  dt  =  - - -  (36) 

0  (-a+iw)^  +  Wq 


The  total  energy  scattered  by  and  transmitted  through  the  aperture  is  then  given  by 


1eIc'^^\ 

^0 


/ 


5(0)) 


•■-«(«  ^-0)^+0)^^  +  4fl^o) 


.2,,2 


-dm 


(37) 


The  integral  in  (37)  can  be  evaluated  using  contour  integration  (Appendix  A)  yielding  the 
following  final  result 
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(38) 


pS  2  2 

ItEnC  (On 


5((OQ+ia) 


a 


-  iS'iw^+ia) 


2.4.4  A  Bound  on  Coupling  for  a  Rectangular  RF  Pulse: 

The  penetration  of  rectangular  RF  pulse  through  an  aperture  and  subsequent  coupling  is 
of  considerable  interest  because  of  the  potential  for  threat  from  high  power  microwaves  (HPM) 
and  EM  interference  from  on-board  or  nearby  radars.  A  bound  on  the  coupling  for  a  rectangular 
RF  pulse  has  not  been  presented  previously  in  the  literature. 

A  rectangular  RF  pulse  can  be  represented  in  the  time  domain  by 

E(f)  =  Eq  coscOq^  (39) 

=  0  otherwise 


The  spectrum  of  this  waveform  is  given  by 


£(o))  =  fE(t)  dt  =  -i£o 


sinfl((o-coQ)  smfl((o+cO(j) 


(o-(o„ 


(O+CO 


0 


(40) 


Substitution  of  this  expression  into  (15)  gives 


W. 


lEy 


tot 


no 

CO 

'/ 


Re 


I 


sm^a((j3 -Of) 

((0-0>o)^ 


S(o))  do 


sin^a(o)  +(i)(j) 


S(o)  do 


((0+(0(^^ 
sma((o-&)n)  sina(o)+o)n) 

/ - r-2 - ^ 

O^-Or 


do 


(41) 
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With  some  difficulty  (Appendix  B),  the  integrals  in  the  above  expression  can  be  evaluated 
using  contour  integration.  This  leads  to  the  following  expression  for 


2jtc^£’o 

^0 


sin2a(Dn 

Re 

a5(o)o)  + 

a+  ° 

• 

0)o  J 

4 

(42) 


Noting  that  i?6[5(<o)]  is  an  even  function  allows  to  be  rewritten  as 


Itic^Eq 

^0 


Re 


/  sm2aci)n) 

2a  +  ® 

5((Oo) 

LV  Mo  ) 

J 

(43) 


Substitution  of  equation  (19)  for  a,  into  (40)  then  yields 


W. 


o,((Oo)Eo 


2r 


tot 


2t1c 


2a  + 


sin2a(<>. 


0)„ 


= 


1  + 


sin2aci>. 


2ao> 


0  ■' 


(44) 


2a(d^ 

Figure  3  -  Influence  of  Pulse  Width  on  the  Normalized  Energy 
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Equation  (44)  has  a  clear  physical  interpretation.  The  term  outside  the  first  pair  of 
brackets  in  (44)  is  the  average  power  that  would  be  transmitted  through  the  aperture  by  a  CW 
signal  at  frequency  Wj,  and  having  amplitude  Eq  .  The  first  term  within  the  brackets,  2a,  is  the 
pulse  length  and,  when  multiplied  by  the  average  power,  gives  a  first  order  approximation  to 
the  total  energy  transmitted  through  the  aperture.  The  second  term  represents  a  correction  to 
the  first  term  to  account  for  the  distribution  of  frequencies  in  a  finite  duration  pulse.  The 
influence  of  pulse  width  on  the  normalised  energy,  =  W^JW^,  is  shown  in  Figure  3.  If 
the  pulse  is  long,  the  correction  becomes  negligible. 

2.4.5  A  Bound  on  Coupling  for  an  Inverse  Exponential  (New  NATOl  Pulse 


While  the  double  exponential  pulse(described  above)  is  often  used  to  describe  a  nuclear 
EMP  because  of  its  simplicity,  it  is  discontinuous  at  time  t=0  and  has  a  singularity  in  its  first 
derivative  at  this  time.  This  discontinuity  results  in  an  overestimation  of  the  high  frequency 
content  in  the  pulse.  To  overcome  this  difficulty,  the  following  inverse  double  exponential 
wave-form  has  been  proposed  [19]. 

E(t)  =  E,-^^-—  .  (45) 

1  +  e 

The  parameters  needed  to  match  the  Bell  Laboratory  double  exponential  pulse  are: 

a  =  10.30  X  10*  s'^ 

P  =  10.34  X  10*  s-‘ 
d  =  1.16  X  10-^ 

=  20  ns 
Eq  =  50  kV/m 


The  spectrum  of  this  pulse  is  given  [4]  by 


E(a)  = 


sinj^(a +i{i))  j 


(46) 


The  total  energy  scattered  by  and  transmitted  through  the  aperture  is  then 
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-2a 


oo 


iSL-2 


W  = 

"tot 


4c^Eqt:\P -a)  ^  a  ^ 
^0 


Re(Siti>)) 


-  cos(^) 


-da 


(47) 


cos] 


The  integral  in  (47)  has  been  evaluated  (Appendix  C)  by  expanding  the  scattering 
amplitude,  S(a) ,  as  a  Maclaurin  series  in  o) .  The  integrals  corresponding  to  the  first  three 
terms  in  this  expansion  were  evaluated  using  contour  integration.  At  low  frequencies,  it  is 
known  [16]  that  the  cross-section  varies  as  o'*  indicating  that  the  third  term  in  the  Maclaurin 
expansion  is  the  first  that  should  be  retained.  Substitution  of  an  expression  of  this  form  into  (47) 
gives  results  that  are  equivalent  to  those  obtained  by  Wame  and  Chen  [4].  At  high  frequencies, 
on  the  other  hand,  the  cross-section  approaches  the  geometric  area  of  the  aperture  and  it  is 
independent  of  frequency.  In  this  case  the  first  term  in  the  expansion  will  be  dominant  and  the 
following  approximate  result  is  obtained  for 


4c'^EoT:\p-a)^  a  P  (2a-p) 


2«  2 


w. 


tot 


UnSin 


2a 


-1 
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ReiSm  . 


(48) 


3.  CAVITY  BACKED  APERTURES 


3.1  INTRODUCTION 

In  a  great  many  cases  of  practical  interest,  sensitive  electronics  are  located  within  a  metal 
structure  (cavity)  that  contains  one  or  more  apertures.  Examples  include  the  avionics  inside  the 
fuselage  of  an  aircraft,  electronic  equipment  in  the  bridge  of  a  ship  or  a  circuit  board  inside  a 
metallic  case.  Often  we  are  required  to  know  if  external  fields  incident  on  these  cavities  will 
interfere  with  proper  operation  of  the  enclosed  electronic  systems. 

In  this  section  of  the  report,  we  examine  the  use  of  a  power  conservation  approach  to 
estimate  the  interior  fields  and  to  place  a  bound  on  the  coupling  of  these  fields  to  the  electronics 
inside.  This  discussion  largely  follows  the  recent  work  of  Hill  et  al  [20]  and  Lee  and  Yang  [18]. 

3.2  POWER  CONSERVATION  AND  STEADY  STATE  CONDITIONS 

Consider  a  plane  wave  incident  on  a  structure  (cavity)  with  apertures  as  shown  in  Figure 
4.  Initially  we  assume  that  the  cavity  is  lossless  except  for  energy  reradiated  back  through  the 
aperture.  This  implies  that  the  cavity  walls  are  perfectly  conducting  and  that  there  is  no  other 
absorption  of  energy  within  the  cavity.  In  this  case,  steady  state  fields  are  achieved  when  the 
power,  P^,  being  transmitted  into  the  cavity  is  balanced  by  the  power,  P^,  being  reradiated  out 
of  the  cavity,  that  is  when 

P,  =  P,  .  (49) 

If  the  total  transmission  cross-section  of  the  aperture  is  o,,  then  the  power,  P,, 
transmitted  into  the  cavity  is 

P,  =  (50) 

where  5.  is  the  power  density  of  the  incident  wave.  If  there  are  N  apertures  in  the  wall  of  the 
cavity,  can  be  written  as  the  sum 

o.'Eo,,-  (51) 

i=l 
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Figure  4  -  Excitation  of  a  Cavity  with  an  Aperture 


A  relationship  similar  to  (50)  can  be  written  for  the  power  reradiated  from  the  cavity. 
In  many  cases,  it  is  necessary  to  express  the  fields  inside  the  cavity  as  a  superposition  of  plane 
waves  and  the  cross-section  that  needs  to  be  used  is  o^,  an  average  over  this  superposition.  In 
general,  the  power  leaking  out  of  the  cavity  is  given  by 

Pr  =  o,  (52) 

where  is  the  power  density  of  the  field  incident  on  the  interior  side  of  the  aperture. 

In  certain  limiting  situations,  the  above  power  conservation  relationship  can  be  used  to 
determine  the  fields  in  the  interior  of  the  cavity  and  the  Q-factor  of  the  cavity.  Some  special 
situations  are  examined  below. 
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3.2.1  High  Freouencv  Approximation 


At  high  frequency,  the  field  distribution  in  the  cavity  can  be  considered  approximately 
uniform  and  the  steady-state  energy,  f/  ,  can  be  written  as 

U^  =  WV  (53) 


where  W  is  the  average  energy  density  and  V  is  the  cavity  volume.  The  average  field  inside  the 
cavity  is  related  to  the  energy  density  by 

W  =  .  (54) 

In  (54),  is  the  permittivity  of  free  space  and  the  electric  and  magnetic  energies  are 
assumed  equal.  The  average  power  density  in  the  cavity,  S^,  is  given  by 

S^  =  ~  =  cW.  (55) 

no 

At  high  frequencies,  the  superposition  can  be  considered  as  consisting  of  plane  waves  of 
all  incident  angles  and  polarizations,  but  only  plane  waves  that  propagate  toward  the  aperture 
contribute  [20]  to  the  leakage.  The  leakage  power  is 

P,  -  <^>|  (56) 

where  (o^)  is  the  average  over  all  incident  angles  and  polarizations  and  is  the  average  power 
density  inside  the  cavity. 


Substitution  of  equations  (50)  and  (56)  into  (49)  yields  the  following  relationship  between 
S^,  the  power  density  inside  the  cavity,  and  S.,  the  power  density  of  the  incident  wave. 


(57) 


The  shielding  effectiveness  of  the  cavity  can  be  defined  in  terms  of  the  ratio  of  the 
incident  and  cavity  power  densities  by  the  following  relationship 
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This  shows  that  the  shielding  effectiveness  of  the  cavity  is  determined  solely  by  the  ratio 
of  the  two  aperture  cross-sections.  In  many  cases  is  larger  than  (a^>  and  hence  the  shielding 
effectiveness  is  negative  indicating  that  the  interior  fields  are  greater  than  the  incident  fields. 


The  Q  of  the  cavity  can  be  determined  from  the  above  relationships.  Noting  that  theQ 
of  the  cavity  is  defined  as 


<?  = 


2tz  X  steady  state  energy  in  cavity 
energy  loss  per  cycle 


Pr 


(59) 


where  «  is  the  steady  state  frequency  and  is  the  steady  state  energy,  then  this  relationship, 
together  with  the  steady  state  power  condition,  can  be  used  to  derive  the  following  expression 
for  cavity  Q 

^  ^  ‘ 

This  shows  that  the  average  transmission  cross-section,  (a^) ,  for  reradiation  determines 
the  cavity  Q  and  that  the  cavity  Q  increases  as  the  aperture  size  decreases.  The  above 
discussion  assumes,  however,  that  the  cavity  is  lossless.  In  a  practical  situation,  wall  losses  and 
absorption  by  systems  inside  the  body  will  reduce  cavity  Q  and  need  to  be  taken  into  account 
if  interior  fields  are  going  to  be  determined  accurately.  The  lossless  situation  described  in  this 
section  corresponds  closely  however  to  numerical  models  that  are  normally  used  to  calculate 
interior  fields. 


A  case  of  interest  is  that  of  uniform  random  excitation  as  discussed  by  Hill  et  al  [20] . 
In  this  case  the  averaged  transmission  cross-section  is  equal  to  the  averaged  cross-section  for 
reradiation  and  the  shielding  effectiveness  is  0  dB. 

Expressions  for  the  aperture  cross-sections  have  been  given  in  the  literature  for  the  cases 
where  the  aperture  is  electrically  large  or  electrically  small.  In  the  following  discussion, 
estimates  of  the  shielding  effectiveness  are  given  for  these  two  situations. 

If  the  aperture  is  electrically  large,  the  aperture  cross-section  can  be  written,  from  a 
geometrical  optics  approximation,  as 


where  A  is  the  aperture  area  and  6'  is  the  angle  of  incidence  of  the  plane  wave  with  respect  to 
a  vector  normal  to  the  aperture. 


The  average  cross-section  for  reradiation  can  be  written  as 

JK 

In  2 

(a  }  =  —f  d4>'  fA  cos0‘  sin0‘  J0‘  =  — 
In  J  J  '> 


0  0 


(62) 


Substitution  of  these  expressions  into  (58)  yields  the  following  for  the  shielding 
effectiveness  of  the  cavity 

SE  =  10  logf — —]  .  (63) 

M  cos0*>' 

This  shows  that  at  normal  incidence  (0‘  =  0)  the  shielding  effectiveness  of  the  cavity  is  - 
6dB  whereas  at  an  angle  of  75°  the  shielding  effectiveness  is  0  dB. 

The  transmission  cross-sections  for  an  electrically  small  circular  aperture  of  radius  a 
were  given  in  section  2.3.  Following  Hill  et  al  [20],  the  averaged  reradiation  cross-seetion, 
(o^) ,  can  be  written  as 

jt/2 

(a  )  =  —  [(  o“  +  Of  )  sinO'  dd'  =  —k*a^  .  (64) 

"  2i  ‘  *  9ti 


For  normal  incidence,  the  shielding  effectiveness  becomes 


SE  =  10 


-L2dB  . 


(65) 


The  above  relationships  indicate  that,  at  high  frequencies  (ie.  frequencies  well  above  the 
lowest  resonance  frequency  of  the  cavity  when  the  field  distribution  inside  the  cavity  can  be 
considered  uniform)  the  interior  fields  are  expected  to  be  of  the  same  order  as  the  incident  field 
—  at  most  a  factor  of  two  above  incident.  It  is  also  of  interest  to  note  that,  at  this  level  of 
approximation,  the  fields  are  independent  of  aperture  size  and  hence  of  cavity  Q.  This 
conclusion  is  supported  by  numerical  calculations  and  measurements  made,  at  DREO,  of  various 
structures  including  the  CPF  bridge  [1]  and  the  electronics  box  of  missile  [24]  (Figure  5). 


3.2.2  Interior  Fields  at  Cavity  Resonance 


For  some  simple  geometries,  the  field  structure  inside  the  cavity  at  resonance  can  be 
expressed  analytically.  For  example,  the  eigenvectors  for  the  lowest  TM  mode  of  a  circular 
cylindrical  cavity  are  given  [10]  by 


(66) 


where  and  are  Bessel  functions  of  the  first  kind,  r  is  the  distance  from  the  centre  of  the 
cavity  and  b  is  the  radius  of  the  cavity.  If  we  consider  a  plane  wave  normally  incident  on  the 
endplate  of  a  circular  cylindrical  cavity  that  contains  a  small,  circular  aperture  as  shown  in 
Figure  6,  then,  by  equating  the  transmitted  and  leakage  powers,  the  following  relationship  is 
obtained  for  the  shielding  effectiveness  of  the  cavity 


SE  =  10 


log 


4 


\|J-<,(2.405^)|2  +  4|J,(2.405|)p 
D  b 


(67) 


where  d  is  the  distance  of  the  centre  of  the  aperture  from  the  centre  of  the  cavity.  The  effect 
of  position  of  the  aperture  on  the  shielding  effectiveness  is  illustrated  in  Figure  7. 


Figure  6  -  Penetration  into  a  Circular  Cylindrical  Cavity  with  Aperture 
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Figure  7  -  Effect  of  Shielding  Effectiveness  on  Aperture  Position 


3.2.3  Interior  Fields  at  Low  Frequencies. 

As  discussed  by  Stratton  [21],  the  dominant  fields  in  the  vicinity  of  an  aperture  at  low 
frequency  are  the  inductive  fields  that  decay  as  r  and  r~^ .  These  fields  account  for  the  stored 
energy.  There  is  a  flow  of  energy  into  the  cavity  associated  with  these  terms,  however  this  flow 
averages  to  zero  over  a  complete  cycle.  As  a  first  order  approximation  to  the  fields  inside  the 
cavity  at  low  frequency,  the  fields  from  the  static  magnetic  and  electric  dipole  moments  of  the 
aperture  can  be  used.  This  assumes  that  the  walls  are  far  enough  away  that  they  do  not 
influence  the  results.  Alternatively,  the  more  sophisticated  mode  matching  approaches  of  Senior 
and  Desjardins  [22]  or  Mendez  [23]  can  be  used. 
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3.3  COUPLING  TO  WIRES  INSIDE  A  CAVITY 


Consider  the  situation  shown  in  Figure  8  where  a  plane  wave  is  incident  on  a  cavity  with 
an  aperture  in  its  surface.  Inside  the  cavity  is  a  wire  that  is  connected  to  an  electronic 
component  such  as  an  integrated  circuit  or  transistor.  An  estimate  (bound)  on  the  energy  that 
can  couple  to  the  wire  is  of  interest  to  determine  if  upset  or  damage  is  likely  to  occur. 


Figure  8  -  Coupling  to  a  Wire  Inside  a  Cavity  with  an  Aperture 


Lee  and  Yang  [18]  have  extended  the  power  conservation  approach  outlined  in  section 
3.2  to  address  this  situation.  Assuming  that  the  field  distribution  inside  the  cavity  is  uniform 
(ie.  high  frequency  approximation),  then  the  power  that  can  be  coupled  to  the  wire  is  given  by 

P  =  a  S  =  a  S.  .  (68) 

w  w  c  p  i 
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In  this  expression  is  the  absorption  cross-section  that  relates  the  power  coupled  to  the 
wire  to  the  interior  field  strength  whereas  is  the  cross-section  that  relates  the  power  coupled 
to  the  wire  to  the  incident  field  strength,  includes  the  effects  of  the  aperture. 

The  Q  of  the  cavity  will  be  modified  by  the  absorption  by  the  wire  and  can  be  written 
more  generally  as 


Q  = 


u>Wj. 

■77 


(69) 


where  Wj.  and  Pj.  are  the  total  energy  and  total  power  dissipated  in  the  cavity  (including  the  loss 
in  the  wire,  the  power  reradiated,  and  other  losses  such  as  wall  losses).  For  conservation  of 
power  ,  the  following  condition  applies 


Using  the  above  relationships,  the  power  density  in  the  cavity  can  be  expressed  as 


aXQ 

S  =  — S. 
‘  2tzV  ' 


(71) 


Using  the  power  conservation  relationship,  one  can  also  write 


Vi 


(72) 


and  thus 


(73) 


As  pointed  out  by  Lee  and  Yang  [18],  the  factor  /Op  is  the  fraction  of  the  total 
power,  Pj,,  being  absorbed  by  the  wire.  If,  for  simplicity,  the  only  loss  mechanisms  are 
absorption  by  the  wire  and  reradiation  through  the  aperture,  then  Op  =  and 


%  = 


0  +0 
w  r 


(74) 
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Certain  limiting  situations  become  clear.  Consider,  for  example,  the  case  when  the 
length  of  the  wire  is  long  compared  with  the  aperture  size  and  .  In  this  case  “  o, 

which  indicates  that  essentially  all  of  the  energy  being  transmitted  through  the  aperture  is 
absorbed  by  the  wire.  In  this  case  the  Q  of  the  cavity  becomes 


Both  the  <?-factor  of  the  cavity  and  the  power  density,  ,  are  much  smaller  in  this  case 
than  the  values  that  would  be  obtained  by  assuming  that  the  cavity  was  lossless.  This  illustrates 
the  hazard  of  using  a  two-step  process  (ie.  firstly,  to  calculate  fields  inside  a  lossless  structure 
and  then,  secondly,  to  use  these  fields  to  estimate  coupling).  It  is  clear  that  the  two-step  process 
will  overestimate  the  coupling  and  violate  the  principle  of  energy  conservation.  If  <  o^, ,  on 
the  other  hand,  then  the  two-step  process  gives  the  correct  result  since  the  wire  does  not 
significantly  load  the  cavity  and  reduce  the  internal  field. 

Lee  and  Yang  [18]  have  used  approximate  expressions  for  the  relevant  cross-sections  in 
(74)  to  determine  a  bound  on  the  coupling  to  a  wire  inside  a  cavity  with  a  narrow  slot  and  have 
compared  these  results  with  experiment.  As  shown  in  Figure  9,  this  procedure  provided  good 
agreement  in  the  cases  reported  by  Lee  and  Yang. 


Figure  9  -  Comparison  of  Measured  Data  [18]  with  Bounding  Estimate 
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4.  SUMMARY  AND  CONCLUSIONS 


This  report  reviews  the  use  of  bounds  to  estimate: 

a)  the  amount  of  energy  that  can  be  transmitted  through  an  aperture; 

b)  the  fields  inside  of  a  cavity  with  one  or  more  apertures;  and 

c)  the  coupling  of  energy  to  wires  inside  of  a  cavity  with  apertures. 

Because  it  is  virtually  impossible  to  model  all  of  the  detail  of  any  real-world  system,  the 
approach  examined  in  this  report  attempts  to  define  the  maximum  energy  or  power  that  can 
couple  to  an  electronic  system  in  terms  of  an  upper  bound  that  is  independent  of  most  of  the 
system  details. 

The  first  portion  of  this  report  reviews  the  methodology  that  has  been  developed  by  Lee 
[2]  and  Wame  and  Chen  [3,4]  to  define  a  bound  on  the  energy  that  can  be  transmitted  through 
an  aperture  when  the  incident  field  takes  the  form  of  a  step  function  or  a  double-exponential 
pulse.  Such  a  bound  is  clearly  of  interest  as,  from  the  principle  of  energy  conservation,  this 
bound  represents  an  upper  limit  to  the  energy  that  can  couple  to  electronic  systems  in  the  interior 
of  the  structure  containing  the  aperture.  Ideally,  however,  one  would  hope  to  be  able  to  define 
a  tighter  bound  to  prevent  overhardening  of  the  equipment  if  the  information  is  used  for  this 
purpose.  We  have  extended  the  work  of  Warne  and  Chen  to  define  bounds  for  the  energy  that 
can  be  coupled  from  a  damped-sine  pulse,  a  rectangular  RF  pulse  and  an  inverse-exponential 
pulse.  These  results  have  not  been  reported  previously  in  the  literature. 

In  the  second  part  of  the  report,  we  show  how  a  power  conservation  approach  can  be 
used  to  estimate  the  fields  inside  a  lossless  cavity  with  apertures.  The  subject  is  of  interest  as 
it  offers  an  alternative  to  the  computationally  intensive  numerical  approach  that  is  often  used. 
Our  work  in  this  area  is  an  extension  of  the  ideas  that  have  been  presented  by  Lee  and  Yang  [18] 
and  Hill  et  al  [20].  Using  the  power  conservation  approach  the  following  conclusions  have  been 
reached: 

a)  at  steady  state,  the  shielding  effectiveness  offered  by  the  cavity  is,  to  first  order, 

independent  of  aperture  size  and  cavity  Q ;  and 

b)  the  minimum  shielding  effectiveness  is  generally  of  the  order  of  -6  dB  provided  the 

frequency  of  the  incident  wave  is  above  cavity  resonance. 
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The  final  part  of  the  report  reviews  the  use  of  the  power  conservation  approach  to  define 
a  closer  bound  on  the  energy  that  can  be  coupled  to  a  wire  or  some  other  absorber  inside  a 
cavity  with  an  aperture.  This  is  essentially  a  summary  of  the  recent  work  of  Lee  and  Yang  [18] . 
Lee  and  Yang  have  shown  that,  by  using  approximate  coupling  cross-sections  for  interior 
absorbers  as  well  as  the  apertures,  it  is  possible  to  define  a  much  closer  bound  than  that  based 
on  the  bound  for  energy  transmitted  through  an  aperture. 

In  conclusion,  we  have  found  that  the  use  of  energy  bounds  represents  a  useful  approach 
for  the  estimation  of  interior  fields  and  coupling  to  wires  inside  a  strucmre  with  apertures.  The 
approach  has  potential  for  simplifying  the  hardening  of  electronic  systems. 
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APPENDIX  A 


In  this  appendix  we  evaluate  the  integral  from  Section  2.4.3: 


■/: 

■/: 


5(ca) 


dui 


_ S((0) _ 

[(o-(a>o  +  ia)]^  [(o-(Wo-ia)]^ 


J(i>  . 


Consider  the  contour  integral 

j  f{z)dz  , 
r 


where 


Kz)  = 


_ W _ 


and  the  contour  F  is 


The  function  /  is  then  analytic  inside  the  contour  F  except  for  the  pole  of  order  2  at 
<jiQ  +  ia  .  The  path  integral  on  the  semicircle  -»•  0  as  R  ^  +  oo .  Thus  hy  the  residue 
theorem 
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2n:i 

lim  M 

(o-tOfl  +  ia 

2ni 

+ 

o 

3 

i  (2xfl)2 

TT 

|5((Do  +  ia) 

2a2 

a 

V 

^(0))  I 

[<o-(o)(,-ia)]2j 

_  25((Jo  +  ia)[ 

(2iay  I 

-  i5'(wQ  +  ia)l  . 


Therefore 


fV...  = 


tot 


no 

TC  Eq  (Oq  J5((0o  +  ia) 


no « 


a 
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APPENDIX  B 


In  this  appendix  we  evaluate  the  three  integrals  from  Section  2.4.4: 


sm^[a((i)  -(Ojj)] 


S((i))  do  , 


h  = 


I*  *°°  sm^[a(ci)  + 

J_„  (<^+«o' 


“o)] 


S(o)  do  , 


and 


I3  = 


/: 


2  sin[a((«) -Wq)]  sin[a((i) +(0(^] 

O^-On 


S(o)  do 


B.l 


Using  the  substitution  x  =  o-o^  and  the  identity 


sin^(aa:)  =  He 


l(l_g2icx^ 

2 


one  obtains 


S(x  +  o)q)  dx  . 


Consider  the  contour  integral 

f  fiz)dz  , 
r 


where 


m 


(l-e2‘«)5(z  +  o)o) 
2z^ 
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and  the  contour  F  is 


The  limit 

1  -p'^^a.z 

lini  z/(z)  =  lim - 5(z  +  (Oo) 

z-o  z~o  2z 

=  -ifl5(6)o)  , 


implies  that 

lim  /  f(.z)dz  =  (x)[-ia5(Oo)](-7i) 
-r(r) 

=  -auSCcOg)  . 


Similar  to  Appendix  A, 


The  contour  F  encloses  no  singularities  of  /  ;  so  by  Cauchy’s  theorem 

f  f(z)dz  =  0  . 


taking  r  ^  0  and  R  -*  +  oo  one  obtains 


f  °‘f(x)dx  -  =  0  . 

J  — oo 


Therefore 


B.2 


=  Re  [«^<S(<«>o)] 


Similarly,  the  integral  I2  ,  with  the  substitution  j:  =  a>+o)o  ,  etc.,  can  be  shown  to 
have  the  value 


5(x-(i)(j)  dx  =  Re  [aTcSC-o)^)]  . 


B.3 


Using  the  identity 


2sin[a(&) -Wq)]  sin[a((i) +o)q)]  =  Re 


g2»a«o_g2itfw 


one  obtains 


=  Re 


g2ia«o_  g2iaa 
(1)^  -Ofl 


5((i))  d(i3  >  . 


Consider  the  contour  integral 

/  f(z)dz  , 
r 

where 


f(z)  = 


^2iaOo_^2iaz 


S(Z)  , 


35 


and  the  contour  F  is 


r  =  [-Rr«o  'r]  u  Y-  o  [-<*>0  -r]  o  y+  [«o  +r,R]  u  r(R) 


Similar  to  the  calculation  in  Rl,  the  limit 


2iaOo_.2i<JZ 

lim  (z  +  Oo)/(z)  =  lim  -  S(z) 


Z--«n 


2isin[2a(i)Q] 

~2ti>rt 


implies  that 


lim  ff(z)dz  =  (i) 


2isin[2a(»)Q] 


-2  Cl)-, 


5(-Wo)(-rt) 


7tsin[2a(i)-j] 


On  the  other  small  semicircle,  the  limit 


lim  (z-Wo)/(z) 
z-+o>0 


lim 

z-+<Do 


^2ia<0(,_^2iaz 

Z  +  COfl 


Siz)  =  0  , 


implies  that 


Also, 
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The  contour  F  encloses  no  singularities  of  /  ;  so  by  Cauchy’s  theorem 


Since 


taking  r  -» 


Therefore 


Finally 


/  /(z)&  =  0 


f  ■  *  /  *  /'  *  /  • 

J  J -/?  J  +r  •'-Wft+r 
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APPENDIX  C 


In  this  appendix  we  evaluate  from  Section  2.4.5: 


4  Eq  tz^  (P  -a)  ^ 


IP...  = 


tot 


^0 


—  -2  — 

'—I. 


Re[S(,(o)] 


oosh(^)  -  oos(l|i) 


d(ji 


A  few  substitutions  will  simplify  the  integral. 


Let  r  = 


2aTi 


-  Tz  ,  then  cos(— -  ^  )  =  -cos/-— 

\  p  /  Ip 

^  ^  /  2x10.3x10^  _  j 
V  10.34x10® 


=  -cosr  .  Note  that 


TZ  «  0.99271  . 


Further,  let  x  -  and  S(x)  =  /?e[S(<i))]  =  /?c[S(-^x)]  .  Since  the  integrand  is  an 

P  27t 


even  function  of  x  , 


A  +09  A  +00 
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and  the  Maclaurin  series  of  S  has  only  even  powers  of  x  : 
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The  problem  is  then  reduced  to  the  evaluation  of  integrals  of  the  form 

^2k 


n  +o> 

J  0 


coshjc  +  cost 


dx  , 


k  =  0,1,2,... 


Let  us  use  the  integral 
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as  an  illustration.  Now 
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and  the  contour  F  is 
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The  function  /  has  two  simple  poles  inside  the  contour  F  at  i(Tr-t)  and  i(7t+t)  . 
This  is  because 


+  2cosf  +  1  =  0  ; 

=  -cost  ±  isinf 
=  cos(7tTf)  +  isin(7iTf)  ; 
z  =  i(7r  -t)  ,  i(v:  +t)  . 


Note  that 


and  the  integrals  on  the  vertical  paths  ^0  as  R  -^  +  oo .  Thus 
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The  residues  are 
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Therefore  by  the  residue  theorem 
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Finally, 


lo  =  lim 
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Using  a  similar  technique,  and  taking  extreme  care  of  the  algebra,  one  obtains 
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Taking  the  first  three  terms  of  the  Maclaurin  series  of  S  ,  one  then  has 
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Since  r  =  0.992  7t  ,  the  high-order  terms  are  relatively  small.  Thus  taking  only  the  first  term, 
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